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[ Abstract. We obtain the nuclear proximity potential by using semiclassical extended Thomas 

Fermi (ETF) approach in Skyrme energy density formalism (SEDF), and use it in the extended 
^-summed Wong formula under frozen density approximation. This method has the advantage 
of allowing the use of different Skyrme forces, giving different barriers. Thus, for a given 
reaction, we could choose a Skyrme force with proper barrier characteristics, not-requiring extra 
"barrier lowering" or "barrier narrowing" for a best fit to data. For the ®^Ni+^'"^Mo reaction, 
the ^-summed Wong formula, with effects of deformations and orientations of nuclei included, 
fits the fusion-evaporation cross section data exactly for the force GSkI, requiring additional 
barrier modifications for forces SIII and SV. However, the same for other similar reactions, 
like '^'^^Ni-l-^^'^^Ni, fits the data best for SIII force. Hence, the barrier modification effects in 
^-summed Wong expression depends on the choice of Skyrme force in extended ETF method. 



1. Introduction 

The unexpected behavior of some fusion-evaporation cross sections at energies far below 
I the Coulomb barrier, has challenged the theoretical models to explain the, so called, fusion 

hinderance phenomena in true coupled-channels calculations (ccc) for reactions such as 
5®Ni-h^®Ni, 64Ni+64]v^j^ j^j^^ 64j^i_^iooj^Q Py_ 'pj^g could, however, be sensitive to the so 
far unobserved, hence not-included, high-lying states. Misicu and Esbensen [2] were the first 
who succeeded in describing the above said three reactions in terms of a density-dependent M3Y 
5^ I interaction, modified by adding a repulsive core potential [3j. The repulsive core changes the 

shape of the inner part of the potential in terms of a thicker barrier (reduced curvature hco) 
and shallower pocket. Here, deformations are included up to hexadecapole (/32 — I^a) and the 
orientation degrees of freedom is integrated over all the allowed values in the same plane. 

The dynamical cluster-decay model (DCM) of preformed clusters by Gupta and collaborators 
[U [5] is found recently [6l [7] to have barrier modification effects as the inbuilt property, where 
"barrier lowering" at sub-barrier energies arise in a natural way in its fitting of the only parameter 
of model, the neck-length parameter. The difference of actually-calculated barrier from the 
actually- used barrier height, corresponding to the neck- length parameter for best-fitted fusion- 
evaporation cross section, gives the "barrier lowering" in DCM, whose values are found to 
increase as the incident energy decreases to sub-barrier energies. Calculations are based on (32 
deformations and orientation 0j-dependent nuclear proximity potential of Blocki et al. [8j. 

Very recently, the Wong formula [9] is also extended by Gupta and collaborators [ini E] 
to include its angular momentum ^-summation explicitly, which is also shown to contain the 



barrier modification effects due to tlie ^-dependent barriers. However, for the ' Ni-based 
fusion-evaporation cross sections [Ij, a further modification of barriers is found essential for 
below-barrier energies, which is implemented empirically either by "lowering the barrier" or 
"narrowing the barrier curvature" by a fixed amount for all ^'s in the potential calculated by 
using the proximity potential of Blocki et al. [8j, but with multipole deformations (32 — (^4 and 
0j-integrated for co-planer nuclei. Apparently, the depth of the potential pocket plays no role, in 
both the DCM and ^-summed Wong formula (the two models are same for capture reactions). 

In this contribution, we use within the ^-summed Wong model, the nuclear proximity potential 
obtained recently [T^] for the Skyrme nucleus-nucleus interaction in the semiclassical ETF 
approach. Using SEDF, the universal function of proximity potential is obtained as a sum of 
the parametrized spin-orbit-density-independent and the spin-orbit-density-dependent universal 
functions (UF's), with different parameters of UF's obtained for different Skyrme forces |12j . 
This method has the advantage of introducing the barrier modifications at sub-barrier energies, 
if needed, by either (i) modifying the Fermi density parameters (the half-density radii and/ or 
surface thicknesses, for "exact" SEDF calculations [E]), (ii) the constants of the parametrized 
UF's [12j or (iii) change the Skyrme force itself since a different Skyrme force would give different 
barrier characteristics (height and curvature). This later possibility is exploited here in this 
paper. It is possible that some Skyrme force would fit the data for one reaction, but not for 
another reaction and hence requiring "barrier modification" or another Skyrme force. 

Section 2 gives briefly the semiclassical ETF method using SEDF, including details of 
approximations used for adding densities. Section 3 discusses the ^-summed Wong formula 
[lOj . Our calculations are given in section 4, and a brief summary of results in section 5. 



2. The semiclassical extended Thomas Fermi (ETF) model 

The SEDF in semiclassical ETF method provides a convenient way for calculating the interaction 
potential between two nuclei. In the Hamiltnian density, the kinetic energy density r as well 
as the spin-orbit density J are functions of the nucleon density pq, q = n,p. For the composite 
system, the densities can be added in either adiabatic or sudden approximation, but we are 
interested in sudden densites since the different terms of Skyrme Hamiltonian density are then 
found to constitute the nuclear proximity potential [121 1131 E! ■ The sudden densities are deflned 
with or without exchange effects (due to anti-symmetrization), and the one without exchange 
effects is also refered to as frozen density [T^. In ETF method, the lowest order r is the Thomas 
Fermi (TF) kinetic energy density tj-f, which already contains a large part of the exchange 
effects, and that the higher order terms include exchange effects in full. Here we limit r(r) and 
J(r) to second order terms for reasons of being enough for numerical convergence [16j. 

The nucleus-nucleus interaction potential in SEDF, based on semiclassical ETF model, is 



VNiR) ^ E{R) - £:(oo) = J H{r)dr 
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where the Skyrme Hamiltonian density 



Hip,r,J) . ^r+'-to 



(l + ^-^o)p'-(xo + i)(p^, + P^) 



12 



(l + ^a;3)p'-(x3 + i)(p2 +p2) 



1 

'4 
1 

16 
1 



In / In 

ti[xi + -) -^2(2:2 + 2) 



(PnTn + PpTp) + — 



3h{l + ^Xi) - t2{l + ^X2) 



In . 1 

iti{xi + -) + t2{x2 + - 



(Vp„)2 + (Vp^)2 
pV ■ J + PnV ■ Jn + PpV ■ Jp 



(2) 



Here, p = pn + Pp, t = Tn + Tp, J = Jn + Jp are the nuclear, kinetic energy and spin-orbit 
densities, respectively, m is the nucleon mass, xi, ti (i=0,l,2,3), a and Wq are the Skyrme 
force parameters, fitted by different authors to ground state properties of various nuclei (see, 
e.g., [13 Us])- Of the available forces, we use the old, well known SIII and SV forces. Coulomb 
effects are added directly. Recently, Agrawal et. al. [T9| modified the Hamiltonian density ([2]) 
on two accounts, and obtained a new force GSKI: (i) the third term in ([2]) is replaced as 



3 r ^ 1 
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and (ii) a new term due to tensor coupling with spin and gradient is added as 



j^ihXi + t2X2).P + ^itl- t2)(J^ + Jt)- (4) 



We have also used this new GSKI force, with additional six, two each of X3, and a, constants. 
The kinetic energy density in ETF method, up to second order [16], for q=n or p is 



5 36 3 6 fq 
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with fq as the effective mass form factor, 

„ 2ml r , xi , , a;2^~l 2to 1 ( , 1, , 1,1 

/,(f) = \^—- + ^) + t2(l + y )} p(r') - ^- |ti(a;i + -) - t^ix^ + 2) j P^l^)- (6) 

Note that both and are each functions of pq and/ or p only. 

The spin J is a purely quantal property, and hence has no contribution in the lowest (TF) 
order. However, at the ETF level, the second order contribution gives 



2m 11^ 

•^9(0 = --T9^0-pPgV(p + p,). (7) 
n ^ Jq 



Note, Jq is also a function of pq and/ or p alone. 

Next, for the proximity potential we introduce the slab approximation of semi-infinite nuclear 
matter with surfaces parallel to x — y plane, moving in z-direction, and separated by distance s 
having minimum value sq. Then, following Blocki et al. |8] and Gupta et al. \12\ [T3l [T^ . the 
interaction potential V]\f{R) between two nuclei separated hy R = Ri + R2 + s, is given as 

/oo 
e{s)ds 

= 27rR J {H{p,T,J)-[Hi{pi,n,Ji) + H2{p2,T2,J2)]}dz 

= 47TRjb(j){D). (8) 

where R is the mean curvature radius, and e{s) is the interaction energy per unit area between 
the flat slabs giving the universal function <j){D) in terms of a dimensionless variable D = s/b, 
with sufrace width 6=0.99 fm. The nuclear surface energy constant 7 = 0.9517[1-1.7826(^^)^] 
MeV fm~'^. (^{D) can be calculated "exactly" or parametrized in terms of exponential and/or 



polynomial functions [12j. For axially deformed and oriented nuclei, R is given in terms of the 
radii of curvature Rn and Ri2 in the principal planes of curvature of each of the two nuclei 
(i=l,2) at the points of closest approach (defining sq), by 
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Here, $ is the azimuthal angle between the principal planes of curvature of two nuclei ($=0^ 
for co-planar nuclei). The four principal radii of curvature are given in terms of radii Ri{ai) and 
their first and second order derivatives i?^(ai) and R'-{ai) w.r.t. Oi, where 

R^{ai) = Ro^[l + J2^^XiYx''\o^^)h (10) 

A 

with i?oi as the spherical or half-density nuclear radius, A=2,3,4..., as multipole deformations, 
and Q!j as an angle between radius vector Ri[ai) and symmetry axis, measured clockwise from 
symmetry axis. For the estimation of sq, we refer to pOj for <l>=0 and to [21j for $ 7^0. 

For nuclear denesity pi of each nucleus (z=l,2), we use the temperature T-dependent, two- 
parameter Fermi density (FD) distribution, which for the slab approximation is given by 



Pi{Zi) = PQi{T) 



1 + exp 



R^{T) 



00 < z < 00 (11) 



a^{T) 

with Z2 = R-zi = {Ri{ai)+R2{a2)+s)-zi, and central density poi (T) = 4^^^(r) [1 + ^\t) ■ 
Then, since pi = + pp., following our earlier work [I2J, for nuclen density we define 

Pn, = {Ni/Ai)pi and p^^ = {Zi/Ai)pi, (12) 



with half density radii Rqi and surface thickness parameters aoi in Eq. (jlip at r=0, obtained 
by fitting the experimental data to respective polynomials in nuclear mass region ^4=4-238, as 

i?o^(T = 0) = 0.9543 + 0.0994^i - 9.8851 x IQ-'^vl^ + 4.8399 x 10"^ - 8.4366 x lO^^^f (13) 

aoi(T = 0) = 0.3719 + 0.0086yli - 1.1898 x lO'^^^ + 6.1678 x 10"^^^ - 1.0721 x 10"^^^ (14) 
The T-dependence in the above formulas are then introduced as in Ref. [22], 

i?o*(r) = i?o*(r = o)[i + o.ooo5T2], aoi(r) = ao»(r = o)[i + o.oir2]. (15) 

Also, the surface width h is made T-dependent [23], 6(T) = 0. 99(1 + 0. 009T^), where T is related 
to the incoming center-of-mass energy Ec.m. or the compound nucleus excitation energy E^jq 
via the entrance channel Qi„-value, as -E^Af ~ ^c.m. + Qin = 5^^^ ~T (T in MeV). 

Next, for the composite system, p = pi + P2, and the t(/>) and J(p) are added as per 
prescrption used. For sudden approximation (with exchange effects), 

t{p) ^ T{pin + P2n) +T{pip + P2p), and J (p) ^ J {pin + P2n) + J {pip + P2p) , (16) 

and for the frozen approximation (equivalently, sudden without exchange effects), 

r{p) = Tiipi) +T2{p2), and J(/>) = /i(pi) + ^2(^2), (17) 

with Pi = Pin + Pip, Ti{pi) = Tin{pin) + Tip{pip), and Ji{pi) = JiniPin) + Jip{pip)- In the following, 
we consider only frozen densities (sudden without exchange), since ^-summed Wong formula 
within sudden approximation does not fit the ®^Ni-based reactions data for use of SIII force |2.4J. 



Finally, adding the Coulomb and centrifugal interactions to the nuclear interaction potential 
Vn{R), we get the total interaction potential for deformed and oriented nuclei [HES], as 

Ve{R) = Vn{R, Ai,|3x^,T, Oi, ^) + Vc{R, Zi,l3x.,,T, 9^,, + Ve{R, Zi,Px^,T, 9^, (18) 

with non-sticking moment-of-inertia Ins {=l^R^) for V^. Eq. (|18|) gives the barrier height Vg, 
position i?^, and the curvature huii for each £, to be used in extended -^-summed Wong's formula 
[lOj . discussed in the following section. 

3. Extended ^-summed Wong Formula 

According to Wong [9j, in terms of ^ partial waves, the fusion cross-section for two deformed 
and oriented nuclei colliding with Ec.m. is 

aiE,,m.,9i,^) = ^J2{2i + l)PeiE,,^.,9i,^), k = ^ -^^^^ (19) 

with ^ as the reduced mass. Here, Pi is the transmission coefficient for each i which describes 
the penetration of barrier Vi{R, Ec.m.,di,^)- Using Hill- Wheeler [26] approximation, the 
penetrability P^, in terms of its barrier height V^{Ec.m., di, ^) and curvature huj£{Ec.m.:Si: ^)) is 

-1 

(20) 

with hiOi evaluated at the barrier position R = i?^ corresponding to Vg. Note, the ^-dependent 
potentials are required here, given by Eq. (jl8p . Carrying out the ^-summation in Eq. (jl9p 
empirically for a best fit to measured cross-section [l.Oj, and on integrating over the angles 
9i and we get the fusion cross-section a{Ec.m)- 

4. Calculations 

We have made our calculations for the ^^Ni+^'^^'Mo reaction, using SHI, SV and GSKI forces, 
with frozen densities. Fig. 1(a) shows for one Ec.m. and fixed (^j, $), a comparison of interaction 
potentials for the three forces, illustrating their barrier height and position to be force-dependent. 
This is an interesting property, which we use to fit the fusion-evaporation cross-section in Fig. 
1(b). Interesting enough, the data [1] for the above said ^^Ni-|-^'"^Mo reaction fit the 0j-integrated 
(<I>=0^) ^-summed Wong formula for only the new force GSKI. The other forces (SHI and SV) 
would apparently need additional barrier modification effects to be added empirically [10]. From 
the deduced ^max-values, presented in Fig. 1(c), we notice that £ 

max B function of Ec.m. vary 

smoothly only for GSKI force, achieving zero value at sub-barrier energies and a tendency to 
saturate at an above-barrier energy. For the other forces (SHI and SV), non-fitting the data, 
the 

^max varies with Ec.m. erratically, in particular at sub-barrier energies, but could also be 
smoothed by adding futher "barrier lowering" or "barrier narrowing" emprically |10j . 

5. Summary and discussion 

Concluding, the ^-summed Wong expression, using the barriers calculated in frozen-density 
approximation in semiclassical extended Thomas Fermi method, based on Skyrme energy density 
formalism, describes the fusion-evaporation cross-section data for ^^Ni-|-^''''Mo reaction nicely 
with the new Syrme force GSKI only. The variation of deduced imax with Ec.m. is found smooth. 
Other Skyrme forces (SHI and SV) demand additional barrier modifications at sub-barrier 
energies for this reaction. However, the same calculation when applied pi] to ^^'^^Ni-l-^^'^^Ni 
data result in a similar good fit, with a smooth dependence of Imax on Ec.m.-, for only the Skyrme 
force SHI. Thus, barrier-modification or no barrier-modification in ^-summed Wong expression 
depends on the choice of Skyrme force in semiclassical ETF method for frozen densities. 
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Figure 1. (a) Interaction potentials for one Ec.m. and fixed {9i, (b) fusion-evaporation 
cross-section as a function of Ec.m. for ^^Ni+^'^'^Mo reaction, calculated by using ^-summed 
Wong formula integrated over 9i {^=0^) and compared with experimental data [1], and (c) 
deduced (max values vs. Ec.m., for the Skyrme forces SIII, SV and GSKI, using frozen densities. 
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